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Interactions and Instabilities in Cosmology's Dark Sector 

Mark Trodden 

Center for Particle Cosmology, Department of Physics and Astronomy, University of Pennsylvania, 
Philadelphia PA 19104, USA. 



I consider couplings between the dark energy and dark matter sectors. I describe how the existence of an adia- 
batic regime, in which the dark energy field instantaneously tracks the minimum of its effective potential, opens 
the door for a catastrophic instability. This adiabatic instability tightly constrains a wide class of interacting 
dark sector models. This talk was presented at, and will appear in the proceedings of the DPF-2009 conference. 



1. Introduction 

Modern cosmology demands the existence of two 
new components to the energy budget of the uni- 
verse. Dark matter is necessary for structure forma- 
tion to occur properly, and to account for numerous 
observations, such as gravitational lensing, the comic 
microwave background (CMB) and galaxy rotation 
curves. The discovery of cosmic acceleration requires 
a second source of new physics, which may come in the 
form of a modification to general relativity, but which 
is perfectly consistent with a cosmological constant 
or a dynamical dark energy component (for reviews 
see flilll). 

At a phcnomcnological level, this description is a 
remarkable fit to all current observations. However, 
at the level of fundamental physics the existence of 
dark matter and dark energy, comprising the vast ma- 
jority of the contents of the universe, poses a critical 
challenge: how do these components fit into our mi- 
crophysical theories of matter and energy? One way 
to search for such relationships is to explore observ- 
able consequences of interactions between dark com- 
ponents and baryonic matter. Examples of this are 
the search for the annihilation products of dark mat- 
ter and collider searches for dark matter. 

Alternatively, if dark matter and dark energy are 
to fit into a unified description, we might expect [|| 
that there would also be interactions between them. 
In this talk I briefly explored some consequences of 
such interactions, focusing on a possible catastrophic 
instability - the adiabatic instability - and a number of 
constraints on stable models. As with all proceedings 
with page limits, thorough referencing is impossible, 
and I have therefore referenced mostly review articles, 
and the papers I actually referred to in the talk. 



where g^ is the Einstein frame metric, is a scalar 
field which acts as dark energy, and are the matter 
fields. The signature is (—,+,+,+) and I define the 
reduced Planck mass by m 2 , = (87rG) _1 . The func- 
tions Oj(</>) are couplings to the j th matter sector. 
The field equations following from this are 



rripG^u 



1 



+ E e4Qj(0) (?) 



3pj) > (3) 



where I have treated the matter field(s) in the ) th sec- 
tor as a fluid with density p 3 and pressure p\ as mea- 
sured in the frame e 2a 'g^ u , and with 4- velocity itj M 
normalized according to g^ u u s ^u sv = — 1. 

For simplicity, let us neglect baryons, and consider 
a composite dark matter sector, with one coupled 
species with density p c and coupling a c (4>) = a(4>), 
and another uncoupled species with density p co and 
coupling a co = 0. Since these both represent types of 
dark matter p c = p co — 0, and it is convenient to set 
Pj = e 3a ' pj to yield 

m 2 p G^ = V M ^V„^- \g»v{V<t>f -V{<t>) gflv 

<4>) 



(4) 



and V M V 



V^ s {4>) = 0. Here the effective poten 
tial is given by V cS {4>) = V(<j>) + e aW p c , and th< 
fluid obeys V Al (p c w^) = 0, and <V„< = ~{g^ 4 
uX)V„a 



2. Modeling Dark Couplings 



I'll focus on the following action, encapsulating 
many models studied in the literature [7[, 

'1 , 1 



S[g ab ,<l>,^j] 



:(V0) 5 



v(4>) 
(i) 



3. The Adiabatic Regime 

The effective potential V c s {<!>) may have a minimum 
resulting from the competition between the two dis- 
tinct terms. If the timescale or lengthscale for <f> to 
adjust to the changing position of the minimum of 
I4ff is shorter than that over which the background 
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Figure 1: The effective potential felt by the dark energy 
field has contributions from the raw potential V(<j>) and 
from the nontrivial coupling ct((f>) of (f> to the background 
dark matter density p. The resulting minimum allows for 
the possibility of an adiabatic regime. 



density changes, the field cf> will adiabatically track 
this minimum |8j. 

In this case the coupled CDM component with cj) 
together acts as a single fluid with an effective energy 
density 



and effective pressure 



Pes(Pc) = ~V[(j) m (p c )} 



(5) 



(G) 



Here (f> m {Pc) is the solution of the algebraic equation 
VM=V'(<P)+a'(<p)e a ^ Pc = (7) 

for ij). Eliminating p c between Eqs. ([5|) and (J6]) gives 
the equation of state p g = p c s(p e s)- 

The coupled fluid acts as the source of cosmic ac- 
celeration, and in the adiabatic approximation the ef- 
fective fluid description is valid for the background 
cosmology and for linear and nonlinear perturbations. 
Therefore, the equation of state of perturbations is the 
same as that of the background cosmology, and the 
matter and scalar field evolve as one effective fluid, 
obeying the usual fluid equations of motion with the 
given effective equation of state. 

A necessary condition for the validity of the adi- 
abatic approximation is that the lengthscales or 
timcscalcs £ over which the density p c varies are large 
compared to inverse of the effective mass of the scalar 
field. A more precise condition can be shown to be 0] 



d]nV[<j> m (p c )] 
din p c 



1 



'cir 



C 2 



< 1 



(8) 



This condition justifies dropping terms involving the 
gradient of <fi from the fluid and Einstein equations. 



Many dark energy models admit regimes in which this 
condition is satisfied for the background and for lin- 
earized perturbations over a range of scales. 

In the adiabatic regime, the inferred dark energy 
equation of state parameter (neglecting baryons) is 



-1 



w 



I — (I — e (a -a))dlny 
\ 'doc 



(9) 



with ctQ = a((j)o) the value today. Thus, w is pre- 
cisely — 1 today, and generically satisfies w < — 1 in 
the past @, [§] 



4. The Adiabatic Instability 

The potential V((f>) may be written as a function 
V(a) of the coupling function a(4>) by eliminating </>. 
This gives, from Eqs. © and (J7|), 



Pcff = V + e a p c 



V 



dV/d(f> 
da/d<j) 



da 



The square of the adiabatic sound speed, c 2 
then given by 



_J_ _ rfjQcff _ dp cS /da 
c 2 a dpcfi dp cS /da 



da 2 
dV 
da 



(10) 



P/p is 



(11) 



In the adiabatic regime the effective sound speed rel- 
evant for local perturbations in pressure and density, 
c 2 s {k,a) = 5P(k, a)/5p(k, a), tends towards the adia- 
batic sound speed and is always negative, since dV/ da 
must be negative so that Eq. ([7]) admits a solution, 
and d 2 V/da 2 must be positive to ensure a positive 

Consider the regime in which this adiabatic limit 
has been reached (c 2 = c 2 ) and focus on a pertur- 
bation with lengthscale C. In order to be in the 
adiabatic regime we require C S> The nega- 

tive sound speed squared will cause an exponential 
growth of the mode, as long as the growth timescale 
~ C/ yj\c 2 \ is short compared to the local gravitational 

timescale m p / y/ p c a(p c )- Since c 2 ml s = (a') 2 V a = 
(a') 2 p e ff/(V/V. a — 1), the instability will operate in 
the range of lengthscales given by 



1 



m c f[(p c ) 



m p \a'[cj> m {p c 
m e s(p c ) 



(12) 



Here the quantity d\nV / da{a) on the right hand side 
is expressed as a function of <f> using a — a((f>), and 
then as a function of the density using <fi — 4> m (p c ). 
In order for this range of scales to be non empty, the 
dimensionless coupling m p |a'| must be large compared 
to unity, i.e., the scalar mediated interaction between 
the dark matter particles must be strong compared to 
gravity (see also [l(| [H| • 
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4.1. Understanding the Instability 

There are two different ways of describing and un- 
derstanding the instability, depending on whether one 
thinks of the scalar-held mediated forces as "gravita- 
tional" or "pressure" forces. 

In the Einstein frame, the instability is independent 
of gravity, since it is present even when the metric 
perturbation due to the fluid can be neglected. In the 
adiabatic regime the acceleration due the scalar held 
is a gradient of a local function of the density, which 
can be thought of as a pressure. The net effect of 
the scalar interaction is to give a contribution to the 
specific enthalpy h(p c ) — J dp/p c of any fluid which 
is independent of the composition of the fluid. If the 
net sound speed squared of the fluid is negative, then 
there exists an instability in accord with our usual 
hydrodynamic intuition. 

On the other hand, in the Jordan frame description, 
the instability most certainly involves gravity. The 
effective Newton's constant describing the interaction 
of dark matter with itself is 



G a r — G 



2m 2 a'(</>) 2 



(13) 



where k is a spatial wavevector |9j. At long length- 
scales the scalar interaction is suppressed and G cc ~ 
G. At short lengthscales, the scalar field is effectively 
massless and G cc asymptotes to a constant. However, 
when m p |a'| 1 there is an intermediate range 



m e s{rn v \a\) 1 < k <C m e ff 



(14) 



over which the effective Newton's constant increases 
like G cc oc k 2 . This interaction behaves just like a 
(negative) pressure in the hydrodynamic equations. 
This explains why the the effect of the scalar interac- 
tion can be thought of as either pressure or gravity in 
the range of scales (fl4| . Note that the range of scales 
(Tj"4")) coincides with with the range (fT2")) derived above, 
up to a logarithmic correction factor. 

From this second, Jordan-frame point of view, the 
instability is simply a Jeans instability. In a cosmo- 
logical background the CDM fractional density per- 
turbation traditionally exhibits power-law growth on 
subhorizon scales because Hubble damping competes 
with the exponential (Jeans) instability one might ex- 
pect on a timescale of 1/^/Gp. In our case, however, 
the gravitational self-interaction of the mode is gov- 
erned by G cc {k) instead of G, and consequently in the 
range (fT4"]) where G cc 3> G the timescale for the Jeans 
instability is much shorter than the Hubble damping 
time. Therefore the Hubble damping is ineffective and 
the Jeans instability causes approximately exponen- 
tial growth. 



5. Examples of Theories with an 
Adiabatic Instability 

1. Exponential Potential and Constant Coupling 

The canonical example is a theory with an expo- 
nential potential of the form 

V = V e- x * /m * , (15) 

with A > and with linear coupling functions 

a{<f) = , (16) 

m p 

where f3 = and C is a constant. The effective 

potential is 



b,p) = V e-^/ m f 



e -/3C<t>/m Pp 



(17) 



and solving for the local minimum of this potential 
yields the relation between and p in the adiabatic 
regime: 



,(A-/3C)<Mp)/m p 



AVb 

-f3C P 



(18) 



Note that C must be negative in order for the effective 
potential to have a local minimum and for an adiabatic 
regime to exist. Restricting attention to this case, 
and defining the dimensionless positive parameter 7 = 
— X//3C, the corresponding effective mass parameter is 



m 



2 \2 -2 T/ 1 + 7 ( P 



(19) 



Using Ijlip we obtain the sound speed squared as 

1 



c? = - 



1+7 ' 



(20) 



so this model is always unstable in the adiabatic 
regime. Eqs. (fl5j) and lfl8]) also yield 



7 



dlnV 



dlnp 1 + 7 



(21) 



which allows one to calculate the range of spatial 
(p) < C < C max (p) over which the in- 
stability operates for a given density p, where 



£ m in(p) 2 



r rag fjVo 

A 2 (l + 7 ) 2 V { p 



^ (22) 



and 



£max(p) 2 = /3 2 C 2 C min (pf 



(23) 



Thus, there is a nonempty unstable regime only when 
j3\C\ 3> 1, ie with the scalar coupling is strong com- 
pared to the gravitational coupling, as we saw earlier. 
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To see the effect of the instability more explicitly, 
consider cosmological perturbations. The Einstein- 
frame FRW equation in the adiabatic limit is 



3m 2 p H 2 = V 



(24) 



where p oc 1/a 3 . This yields a(t) oc t 2 /^ +3wett \ where 
the effective equation of state parameter is 



W c ff 



1 + 7 



In the strong coupling limit |C| — ► oo, w c ff — ■> — 1. 
Thus the adiabatic regime of this model with large 
\C\ is incompatible with observations in the matter 
dominated era, where w c s ~ except for at small 
redshifts. Nevertheless, the model is still useful as an 
illustration of the instability. 

From ([22]) , (j23j) and (|24|) the range of unstable scales 
is given by 



1 



/3 2 C 2 



H 2 a 2 



1 



3(1 + 7) 



(26) 



where k is comoving wavenumber. This range of scales 
always lies just inside the horizon. A given mode k will 
evolve through this unstable region before it exits the 
horizon. 

One can then show that the perturbation evolution 
equation, specialized to the exponential model, and in 
the strong coupling limit \C\ — > oo is 



da 2 



3 dd 
a da 



k 2 



7 S = 



(27) 



In the strong coupling limit H is approximately a con- 
stant, H « Hq, and the growing mode solution is 



6(a) « ~K X 

a \rloa 



2ka 



exp 



H a 



(28) 



where K\ is the modified Bessel function. The mode 
grows by a factor ~ e when the scale factor changes 
from a to a + Aa, where Aa/a ~ aHo/k <C 1 for 
subhorizon modes. 

A more detailed analysis of the cosmology of this 
model is given in [l2l ] , but in the non-adiabatic regime 
\C\ ~ 1 rather than the strong coupling regime \C\ 3> 
1 considered here. 

2. Two Component Dark Matter Models 
As a more realistic example, consider models in which 
there are two dark matter sectors, a density p c which 
is not coupled to the scalar field, and a density p co 
which is coupled with coupling function (|16[) and ex- 
ponential potential ([15]) . Both of these components 
are treated as pressureless fluids. The FRW equation 
for this model in the adiabatic limit is 



3mlH 2 = V 



e p c 



(29) 



The first two terms on the right hand side of Eq. (|29l) 
act like a fluid with equation of state parameter given 
by ([25]) . and in the strong coupling limit |C| » 1 
this fluid acts like a cosmological constant. Thus, the 
background cosmology can be made close to ACDM 
by taking |C| to be large. 

The fraction of dark matter which is coupled must 
be small in the limit of large coupling, \C\ ^> 1. De- 
noting Q v = V/(3mlH 2 ), Q co = e a p co l\Zm\H 2 ) and 



(25) n c = p c /(3m 2 p H 2 ) gives 1 = fi x 



fL. Also 



from Eq. (I18p it follows that, if the asymptotic adia- 
batic regime has been reached, f2 co = "/£lv, yielding 



a. 



7 



1 + 7 



(30) 



Since il c ~ 0.3 today, and 7 <^ 1 in the strong cou- 
pling limit, we must have f2 co <C 1 today. 

The maximum and minimum lengthscales for the 
instability are still given by Eqs. (|22[) and (p3|) . but 
with p replaced by p co . Since p co is approximately 
a constant in the strong coupling limit, these length- 
scales are also constants. If the parameters of the 
model are chosen so that f2 c ~ 1 today, then 



0\C\ 



(31) 



The evolution equations for the fractional density per- 
turbations (5j = 5pj/pj in the adiabatic limit on sub- 
horizon scales are given by 



S c + 2HS C 
S co + 2H5 co 



1 



2m^ 



2m^ 



PcSc 



PcSc 



1 



2m^ 



Ct r 
e PcoOc 



(32) 



2ml 



2(3 2 C 2 



k 2 



e a p co Sc 



(33) 



The condition for the instability to operate is that the 
timescale associated with the second term on the right 
hand side of Eq. (|33p be short compared with if -1 , 



or 



P 2 C 2 k 2 
2 ? P c 



7 e a > E 2 m\ 



(34) 



Now the effective mass for this model is given by 



'off 



(3 2 C 2 m- 2 {V + e a p co ) = 3l3 2 C 2 m~ 2 H 2 (n v 



fl co )- Substituting this into Eq. (f3~4"|) and using Eq. 
([30]) gives the criterion k/(aH) ^> 1. Therefore the 
instability should operate whenever modes are inside 
the horizon and in the range of scales (f3T|) . 

These expectations are confirmed by numerical inte- 
grations. Figure[2]shows a numerical analysis of such a 
two component model with an exponential potential 
with A = 2, strong coupling (C = —20), and typi- 
cal cosmological parameters, Hq = 70 kms -1 Mpc -1 , 
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Figure 2: [Bottom] The two component coupled dark en- 
ergy (CDE) model. At late times the scalar field finds the 
adiabatic minimum with asymptotic equation of state, and 
sound speed = —1/(1 + 7) = —0.89, able to reproduce a 
viable background evolution consistent with supernovae, 
CMB angular diameter distance and BBN expansion his- 
tory constraints. The figure shows the evolution of the 
effective equation of state, w e ff (black full line), the adi- 
abatic speed of sound, c\ = P/p for all components (blue 
long dashed line) and for the coupled components only 
(green dot long dashed line), and effective speed of sound 
for = SP/Sp at k — 0.01/Mpc for all components (red 
dot-dashed line) and for the coupled components alone 
(magenta dotted line). The effective equation of state for 
a comparable ACDM model with fi c = 0.25, fit = 0.05 
and J1a = 0.7 is also shown (black dashed line). [Top] 
The growth of the fractional over-density 8 — Sp/p for 
k = 0.01/Mpc for the coupled CDM component, S co , (red 
long dashed line) and uncoupled component, S c , (black full 
line) in comparison to the growth for the ACDM model 
(black dashed line). 



n b = 0.05, n c = 0.2, fl co = 0.05, and fl v = 0.7. The 
existence of a dynamical attractor renders the final 
evolution largely insensitive to the initial conditions 
for <f>, and the effect of the coupled CDM component 
peculiar velocity in the initial conditions can be ne- 
glected, since it is many orders of magnitude smaller 
than the density perturbation. The background evolu- 
tion is entirely consistent with a ACDM like scenario. 
The large coupling drives the evolution to an adia- 
batic regime at late times, with an adiabatic sound 
speed c\ — > —1/(1 + 7) as in (2D]). This drives a rapid 
growth in over-densities once in the adiabatic regime 
so that although consistent with structure observa- 



tions at early times, they are inconsistent once the 
accelerative regime has begun. 

In summary, these models provide a class of theo- 
ries for which the background cosmology is compati- 
ble with observations, but which are ruled out by the 
adiabatic instability of the perturbations. 



6. Conclusions 

I have briefly discussed a broad class of models in 
which dark matter is coupled to a dark energy compo- 
nent, assumed to be responsible for the acceleration of 
the universe. Within these models, I have discussed 
a possible instability - the adiabatic instability - that 
arises in a range of cosmological and astrophysical set- 
tings, and which rules out a set of parameter values. 
I have discussed specific examples of models in which 
this instability is active, and it is worth noting that 
one may carry out similar analyses to constrain sub- 
classes of chameleon and MaVan models. 
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